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Abstract 

A SIC consists of N 2 equiangular unit vectors in an N dimensional Hilbert 
space. The frame potential is a function of N 2 unit vectors. It has a unique 
global minimum if the vectors form a SIC, and this property has been made 
use of in numerical searches for SICs. When the vectors form an orbit of 
the Heisenberg group the frame potential becomes a function of a single 
fiducial vector. We analytically compute the average of this function over 
Hilbert space. We also compute averages when the fiducial vector is placed 
in certain special subspaces defined by the Clifford group. 
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1. Introduction 



Symmetric Informationally Complete Positive Operator Valued Measures, or 
SICs, is the unwieldy name for a simple idea [1,2]: a set of N 2 unit vectors 
in an N dimensional Hilbert space, equiangular in the sense that 



If the vectors are reinterpreted as projectors, that is as points in the set of 
Hermitian matrices of unit trace (the space where the density matrices live), 
they form a regular simplex in an N 2 — 1 dimensional Euclidean space. This 
also explains why we want their number to be N 2 . The corners of such a 
simplex can be used to define barycentric coordinates for any density matrix, 
which is what "informationally complete" stands for. In quantum informa- 
tion theory SICs have attracted attention because they — if they exist — are 
useful for quantum state tomography [2, 3]. In quantum foundations they 
have attracted attention as a preferred kind of measurement [4] — and they 
have been studied in many other branches of science under names such as 
"equiangular lines" [5], "equiangular tight frames" [6], and "maximal quan- 
tum designs" [1]. Strohmer and Heath provide a mathematical survey [7]. 

The question whether SICs exist for all N has turned out to be extraor- 
dinarily difficult to answer. They have been constructed in most (but not 
all) dimensions N < 19 [1, 8, 9, 10, 11]. Numerical searches have been suc- 
cessful for all N < 45 [2], but no general formula has emerged. This is a 
bit surprising, given that we are really asking a very simple question about 
the shape of the convex body of density matrices: is it possible to inscribe a 
regular simplex in this body, with N 2 corners on its outsphere? The avail- 
able evidence does however suggest that the answer is "yes" for all N, and 
moreover one can always find a SIC covariant under the Heisenberg-Weyl 
group, meaning that it can be constructed by first choosing a fiducial vector 
|^o)> and then acting on this vector with the N 2 elements of the (projective) 
Heisenberg-Weyl group. 

A natural measure of how close a given set of N 2 unit vectors \ipi) is to 
forming a SIC is given by the function 
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By construction / = if and only if the vectors form a SIC. The function / 
is related to the frame potentials 
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These frame potentials also assume their global minimum for a SIC, and in 
fact it is enough if F 2 assumes its minimum, since i*\ is known to follow 
suit. Therefore numerical searches for SICs have focused on minimizing F 2 . 
Actually a frame potential is defined for any integer t [2, 6], but this does 
not concern us here. 

We will have nothing to say about the existence problem here, rather we 
will compute averages of the function /, with and without the assumption of 
group covariance. We use the Fubini-Study measure on complex projective 
space to perform the averaging. We also compute averages when the fidu- 
cial vector is restricted to lie in certain subspaces defined by the Heisenberg 
group and its normalizer (the so-called Clifford group). These calculations 
were made for a reason: in the course of an investigation of certain config- 
urations of N 2 vectors that occur in a related problem (to be precise but 
perhaps not informative, the torsion points of an elliptic curve defined by a 
complete set of mutually unbiased bases) we computed the values of / for 
these configurations. And the question then arose whether the values we ob- 
tained were high or low. To answer questions like this — and there are many 
such questions — one needs to know these averages. 

Our final results have a certain elegance. Our paper is organized as 
follows: in section 2 we state some definitions, and in section 3 we describe the 
special subspaces we are interested in. In section 4 we give details concerning 
our calculations, which were carried out by brute force. Readers who are 
familiar with the SIC problem, and who do not want the details of an involved 
calculation, are advised to go directly to section 5, where we state our results 
and make some comments. 3 

3 The first version of this paper contained a mistake, in eq. (34). A counterintuitive 
and exciting result ensued. The mistake was spotted by Christopher Fuchs and by an 
anonymous referee; it is better to be right than to be exciting, and we thank them for this. 
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2. Definitions 

The Heisenberg-Weyl group [12] is the group generated by two elements r 
and a subject to the relations 

ar = qra, r N = a N = 1 , q = e 2m ' N . (5) 

Up to unitary equivalence there is a unique unitary representation of this 
group such that 

a\a) = q a \a) r\a) = \a + l), 0<a<N-l. (6) 

This is how we fix coordinates in Hilbert space; vectors can then be defined 
by their components 

N-l 

m = e z» ■ (7) 

a=0 

Up to phases, a general group element is 

Aj=tV, 0<i,j<N-l. (8) 

The phases do not matter to us. 

Now consider N 2 unit vectors forming an orbit under the Heisenberg 
group, 

= A,#o> , (9) 

where j^o) is some fiducial unit vector. The frame potential evaluated on 
such an orbit becomes a function on the projective Hilbert space CP^ -1 . 
This is the case we are most interested in, so we define 

N 2 / 1 \ 2 

1* = ^- e K^>r-T— r . (io) 

Because of group covariance there are only iV 2 — 1 terms in the sum. 

The definition of /# can be manipulated further [4, 13]. Following Ap- 
pleby, Dang and Fuchs we use an explicit matrix representation of A? to 
write 
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{M A,#o) = Z a q b ^i +l Z b = q^Z a Z a ^ . (11) 
Then we take the Fourier transform 

I]r^|(V>o|A,#o)| 2 = Y. 2 aZ a+ k-iZ a+ kZ a -i . (12) 
3 a 

From this it is easy to show that the first frame potential Fi = N 3 (this is 
true for all SIC-POVMs), and moreover that 
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(13) 



This is the expression we will work with in the sequel. 

We will be interested in averages of / and To compute these averages 
we use the Fubini-Study measure d//ps on projective Hilbert space [14]; this is 
the natural definition of an average in the absence of any special information. 
What we wish to compute is 

(/> B vol[CP"-'] / d " FS 1 • (14) 

and similarly for f H . Computing (/) is straightforward: using first the linear- 
ity of the expectation value and then the unitary invariance of the measure 
we obtain 

(/)^{(,(«)P-^) ! ). a, 

Here j^o) is any fixed vector. To compute (/#) requires more work — we will 
fall back on the explicit expression (13), and then collect terms. 

We have now defined the functions we wish to average, and we have 
defined "average" . It remains to define the special subspaces of Hilbert space 
that we are about to consider. 



3. Special subspaces 
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We will compute averages of fu also when the fiducial vector is confined to 
lie in certain interesting subspaces of the Hilbert space, picked out by the 
Clifford group. By definition the latter is the normalizer of the Heisenberg- 
Weyl group in the group of all unitaries. Thus the Clifford group is the group 
of all unitaries U such that 

UDijrf = uDi.j, , (16) 

where u is a phase factor. We are interested in representations up to a phase, 
and it can be shown that the Clifford group with the phases ignored is iso- 
morphic to a semi-direct product of a symplectic group with the Heisenberg 
group itself. The importance of this automorphism group was stressed by 
Zauner [1] and Grassl [8]; for a self-contained account see Appleby [9]. 

When the dimension is odd the Clifford group contains elements of order 
2. They play a major role in the definition of discrete Wigner functions 
for the odd dimensional case [15, 16], and for the special case when the 
dimension is an odd prime number they are also known as Wootters' phase 
point operators [17]. They are symmetries of an elliptic curve associated to a 
complete set of mutually unbiased bases [18]. An elliptic curve can be defined 
as an embedding of a torus into complex projective space. The Heisenberg 
group acts on this torus, and gives rise to N 2 torsion points on the curve, 
roughly analogous to the iVth roots of unity on a circle. Hughston [19] has 
made the interesting comment that, in the special case N = 3, these torsion 
points define a SIC. In general it is known that the torsion points lie in the N 2 
subspaces defined by elements of the Clifford group of order 2. In particular 
there is such an element A obeying 

ADijA = D-i-j . (17) 

Because A 2 = 1, the operator A is both unitary and Hermitean. It splits 
Hilbert space into two subspaces Ti^ of dimensions n and n — 1, respectively, 
where N — 2n — 1. Explicitly these subspaces are defined by 
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Unfortunately the fiducial vectors of the SICs do not lie in these subspaces 
when N > 3. However, we believe that /# averaged using the Fubini-Study 
measure in such subspaces gives some feeling for the systematics of the SIC 
problem, and we will compute this average. 

For any N the Clifford group contains elements of order 3. Zauner [1] 
conjectured that there always exists a SIC such that the fiducial vector is an 
eigenvector of a symplectic transformation of order 3, and his conjecture has 
been verified by Appleby [9] in all available cases [2] . We would therefore like 
to know the average of fn over such subspaces. Unfortunately it is not so 
easy to describe these subspaces for arbitrary N, and we therefore confined 
ourselves to the special case N = 7. Then there exists an element of order 3 
acting like a permutation matrix, and the subspaces are explicitly 
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(19) 



where a, the eigenvalue, is a primitive third root of unity. The dimension of 
the subspace Hi is 3, and the dimensions of the two orthogonal subspaces 
are 2. There is a fiducial vector for a SIC in Hi [9]. It seemed reasonable to 
expect that the average f H over this subspace would be quite low, but this 
expectation was not borne out. 



4. Calculations 
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In order to average the frame potential (13) over the entire Hilbert space, 
we observe — see eq. (25) below — that the angular integrals will make all 
non-real terms of the sum disappear. We also note that a term of the form 
\Zi\ A \Z 2 \ A gives the exact same contribution to the average as one of, say, the 
form |Z 3 | 4 |Z 5 | 4 , and all one needs to keep track of is the number of different 
\Zj\ in a term, as well as the exponents present. Consequently, we look at 
the sum (13), and ask in how many ways we can get a term of the type |Zi| 8 , 
in how many ways one of the type |Zi| 6 |Z 2 | 2 , and so on. The problem has 
then been reduced to elementary equation solving. 

Computing the average of fu as given in (13), amounts to computing the 
average 
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(20) 



For calculational convenience, the sum can be split into two parts: the i — 
part where all terms are real and we have to keep track of cross terms, and 
the remaining % ^ part, where cross terms contribute only in special cases: 
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Let us first consider the i — term, which can be rewritten as 

N-l N-l N-l 



(21) 
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k=0 a=0 b=0 



We will divide our analysis into two cases, the case a = b and the case a ^ b. 
When a = b the terms we get will be either of the type \Zi\ 8 — this will occur 
for k = — or of the type \Zi\ 4 \Z2\ A , which will be the case for all remaining 
values of k, no matter what is the value of a and b. The total contribution, 
modulo integration, from the % — term will for a = b then be 



iVlZil 8 + iV(iV — 1)\Z 1 \ 4 \Z 2 



(22) 



where the factor iV in front of the first term comes from the number of 
possible choices of a = b, and the second term is multiplied by the factor 
N(N — 1) for iV different choices of a = b and N — 1 choices of k ^ 0. 
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We now turn to the sum for different values of a and b. In this case, 
we will have four different types of \Zi\, each to the power of two, in each 
term, except in the three special cases where the values of two or more of 
the indices coincide, namely when a = b + k, b = a + k and k — 0. In 
the first two cases we get terms of the \Zi\ A \Z 2 \ 2 \Z^\ 2 type, in the last we 
get a \Zi\ i \Z 2 \ i type term. The remaining choices of k, a and b — there are 
N(N — 1)(N — 3) of them — give terms with all indices different, i.e. terms 
of the type \Zi\ 2 \Z2\ 2 \Z 3 \ 2 \Z i \ 2 . The total contribution from the i — term 
is given in table 1. 

Moving on to the % 7^ case, we can forget about cross terms except when 
k — 0. Consequently, when k 7^ we can consider the sum 

N-l N-l N-l 

X! X! X! \Za\ \Z a +k-i\ \Z a +k\ \Z a -i\ , 
i=l k=l a=0 

while when k = we still have to work with the expression from equation 
(21). As in the previous case, we ask what special cases occur, i.e. when 
two or more indices take on equal values. These cases turn out to be k = i, 
k = —i and k — 0. The respective contributions to the average are given in 
the table, as is the contribution for all other choices of indices. 
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Table 1: Number of different type terms in the average for odd N. 



Subtracting the constant term and multiplying by the initial factor of 
(13) yields an expression which can be integrated over all of Hilbert space 
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in order to obtain the SIC-function average. To perform the integrals we 
parametrize the unit vectors as 

Z a = (Vp~o, y/w?*) , (23) 
where the ranges of the parameters are 

p + Pl + . . . + pn^ = 1 , p > , p k > , < v k < 2tt . (24) 
We solve for p , and obtain the explicit expression [14] 

(_/V— 1)' f 1 rl-pi-...-pN-2 /-27T r2w 

{f) = J^W ZT Jo dPl ---Jo dUl ---Jo dUN ~ lf - 

(25) 

In calculating the averages a number of standard integrals will be used, 
namely 
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(26) 
(27) 
(28) 
(29) 
(30) 

The computations for odd values of N are straightforward, with few com- 
plications. Even N are, however, a somewhat different story. To begin with, 
cross terms play in also for non-zero values of i and k, and in a somewhat 
different manner than in the odd iV case. Whereas we get cross term contri- 
butions for odd N only when all terms are real in a sum — as for k — 0, for 
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instance — we here find pairwise equal terms inside the modulus sign of the 
i 7^ term of equation (21), so that care has to be exercised when squaring. 
Using the same strategy as for odd N — dividing the sum into % — and % ^ 
parts and asking when two or more indices are equal — we find a number of 
index choices for which terms of order higher than two occur, some of which 
coincide with the ones for odd N. As opposed to in the odd N case, though, 
the number of possible choices for the summation indices depend on the par- 
ity of i, as well as on the modulo four value of N, so that the calculational 
details differ beteween, for instance, N = 6 and N = 8. This is because the 
restrictions one gets when solving the equal-indices-equations, are sometimes 
equivalent for, say, % — 3 and N = mod 4 and independent for iV = 2 mod 
4 for the same i, giving a different number of combinatorial possibilities. The 
final results for the frame potential average taken over Hilbert space, on the 
other hand, turns out to be identical for the two cases; this average is given 
in section 5. 

Also when considering the sum average over 7i + and 7i~ , the calcula- 
tions get more involved than in the non-restricted odd N case. As for even 
N, we get cross term contributions to the integral, and further technical 
complications are due to the fact that the respective subspace conditions, 

Zi = Z_i (31) 

and 

Zi = -Z_i , (32) 

make the number of equations to solve when asking what possibilities there 
are for the power of any Zj to be higher than two in a term larger. In other 
words, the calculations have to be divided into a greater number of different 
cases, but the logic is much the same as in the calculations for all of Hilbert 
space sketched above. In using equations (26 - 30), one should keep in mind 
that the dimension of the space over which the average is taken is no longer 
equal to N. 

5. Results 

Our results are easy to state. The function / attains its global minimum at 
a SIC, and its global maximum if all the iV 2 unit vectors coincide. Thus 



11 



TV 4 N - 1 

< / < . (33) 

Its Fubini-Study average value is 

N 2 N -1 

</Vs = TTvTT ' (34) 

When the dimension is large, the scalar product between randomly choosen 
vectors is close to zero, and this is reflected by the average. 

If we specialize to A 2 unit vectors forming an orbit under the Weyl- 
Heisenberg group we believe that the global maximum occurs if the fiducial 
vector is an eigenvector of some element in the group, in which case the orbit 
collapses to an eigenbasis. Thus 



A 3 A 

0<f H < 



2 A + l 



(Unfortunately we were unable to prove that the eigenbasis represents the 
global maximum. It seems obvious that this is so, and we did check that it is 
a local extremum. For safety, we do not number this equation.) The average 
value depends on whether the Hilbert space dimension is odd or even: 

jV 2 N(N-l) 
* = *.-!: U*)- = T {N+ \ KN l l) (35) 

A 2 N 2 

N = 2n : (/„>„ = - (jy - ^ - ^ . (36) 

Asymptotically there is no difference. It is perhaps somewhat unexpected 
that the average of /# has the same asymptotic behaviour as the average of 

We also computed the average when the fiducial vector is restricted to lie 
in certain interesting subspaces defined by elements of the Clifford group (the 
normalizer of the Weyl-Heisenberg group). Elements of order 2 occur in odd 
dimensions N — 2n — 1. Their eigenspaces have dimension dim[7i~] = n — 1 
and dim[7i + ] = n. As observed by Hughston [19], when = 3 H~ contains 
only one vector, and it is a fiducial vector for a SIC. However, 3 is a very 
special dimension. The complementary subspace TC + contains a SIC fiducial 
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vector too, as well as four vectors from a complete set of four mutually 
unbiased bases, for which /# attains its (conjectured) maximum. The average 
over H + is (f H ) = 81/40. For iV > 3 the situation is quite different. We find 

N = 2n - 1 > 3 : (/„)„- = </„> M+ = N* ^'^ ■ (37) 

Asymptotically this is twice the average over the full Hilbert space. It is a 
bit striking that the averages are equal. 

It is harder to deal with the subspaces defined by elements of order 3. 
At the same time they are more interesting, because Zauner's conjecture 
[1, 9, 20] says that the fiducial vector can always be choosen to lie in the 
subspace with eigenvalue 1. We did compute averages for N = 7. There are 
three subspaces Tii, Ti, a , 7i a 2, labelled by the eigenvalues of the element that 
cubes to one (a = e 2m ^). Their dimensions are respectively 3, 2, and 2. The 
average values of fn in the three subspaces are 

151 • 7 3 

(f H )m = 5^33 « 15.985 (38) 
37 • 7 3 

(f H )n a = (f H )n a2 = ^T-ys « n - 751 • ( 39 ) 

The average over the subspace that contains the SIC fiducial vector is above 
the average (fn) over the full Hilbert space. This surprised us. 

The values of fn vary widely, also within the special subspaces (except 
that fn is actually constant on Hr when N — 5). We describe the situation 
for N = 7 in the following table: 





/ 
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fn(n + ) 




him) 


fa (H a ) 


Minimum 








12.2 (?) 


4.764 (?) 





? 


Average 


18.375 


14.29 


25.72 


25.72 


15.98 


11.75 


Maximum 


900.4 


128.6 (?) 


128.6 (?) 


42.88 (?) 


? 


? 



Some entries were left uncomputed, and some without proofs — we marked 
the latter with a question mark within bracket, even though we are sure they 
are correct. 
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Unfortunately we are unable to see how to compute the asymptotic be- 
haviour of (fn) for a fiducial vector in subspaces of the kind considered by 
Zauner. We expect them to depend on number theoretical details of the 
dimension. It would be interesting to see numerical studies of such averages 
however. 
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